








EHEATA NO. 1 


THURTy-SETEHTH ANNUAL REPORT 
OF THE 

NATIONAL ADVISORT COMMITTEE FOR ASIOHAUTICS 

1951 

REPORT 1028 


Page 665, column 1 : In equation (38a), the last bracketed expression 

shoiild he corrected as follows: 

[ip^ + 6 - 4xo(2_ + 32^ _ . - . . 

Page 665, column 2 : In equation (38b), the last bracketed expression 

should he corrected as follows: 


[^2 + ij.)(4 _ 5xQ^ 


Page 665, coltJirm 2 ; In equation (k-^c), the factor 2 preceding the second 
parenthesis should he deleted} that is, the second term, within the 
bracket should read 

- (1 - 2xo)(fi + 

Page 665, column 2 : Equation (i(-6a) should be corrected to read as follows: 


r 


Ml = - S 


it 

pit 


(1 - 2xo)(fi + Gi) + 


+ 2/= . = 


(f3 -h G5) 


+ ^(2 - 3 xq) 


U 3P' 


It is pointed out that the foregoing errors have been corrected in a 
subsequent HACA publication (NACA TN 3O76 by Nelson, Rainey, and 
Watkins ) . 
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EFFECT OF ASPECT RATIO ON THE AIR FORCES AND MOMENTS OF 
HARMONICALLY OSCILLATING THIN RECTANGULAR 
WINGS IN SUPERSONIC POTENTIAL FLOW^ 

By Charles E. Watkikb 


SUMMARY 

I 

This repoti treats the effect of aspect ratio on the air forces | 
and moments of an oscillating flat rectangular wing in super- 
sonic potential flow. The linearized velocity potential for the 
wing undergoing sinusoidal torsional oscillations simultane- 
ously with sinusoidal vertical translations is derived in the 
form of a power series in terms of a frequency parameter. The 
series development is such that the differential equation for the 
velocity potential is satisfied to the required power of the fre- 
quency parameter considered and the linear boundary conditions 
are satisfied exactly. The method of solution can he utilized 
for other plan forms-^that is, plan forms for which certain 
steady-state solutions are known. 

Simple, closed expressions that include the reduced frequency 
to the third power, which is sufficient for application to a large 
class of practicable prohhms, are given for the velocity potential, 
the components of total force and moment coefficients, and the 
components of chordwise section force and moment coefficients. 
The components of total force and moment coefficients indicate 
the over-all effect of aspect ratio on these quantities; however, 
the com ponents of chordwise coefficients yield more information 
because they account for the spanwise distribution of aerody- 
namic loading of a rectangular wing and may therefore be 
more useful for flutter calculations. It is found that the com- 
ponents of force and moment coefficients for a small-aspect- 
ratio wing may deviate considerably from those of an infinite- 
aspect-ratio wing. Thickness effects which may alter some of 
the conclusions are not taken into account in the analysis. 
Results of some selected calculations are presented in several 
figures and discussed. 

INTRODUCTION 

The effect of aspect ratio on the single-degree torsional 
instability of a finite rectangular wing oscillating in a super- 
sonic stream was treated in reference 1 by expanding, in 
powers of the frequency of oscallation, the linearized velocity 
potential developed in reference 2. Since only alow oscilla- 
tions were considered pertinent to single-degree torsional 
instability, terms in the expansion involving the frequency 
of oscillation to powers higher than the first were not 
considered. 

In the present report the expanded linearized velocity 
potential is used to study the effect of aspect ratio on the air 

^ Supersedes XACA TX 2064, “Effect of Aspect Ratio on. tbe Air Fotccs and Moments 


forces and moments of an oscillating, thin, flat, fimte, rec- 
tangular wing when higher powers of the frequency of oscilla- 
tion are taken into account. The motions considered are 
sinusoidal torsional oscillations about a spanwise axis taken 
simultaneously with sinusoidal vertical translations of this _ 
axis. The velocity potential is developed by use of sources . 
and doublets, so as to include aU. powers of the frequency 
of oscillations up to any desired power. Simple, closed 
expressions are given for the velocity potential, components 
of the total force and moment coefficients, and components 

of the chordwise section force and moment coefficients involv- . 

ing powers of the frequency up to and including the third 
power. Extension of the results to include higher powers of 
the frequency is straightforward. 

A recent publication, reference 3, that became available 
after this iuvestigation was completed, is partly devoted to 
the treatment of a rectangular wing undergoing the same 

types of harmonic motions as those considered herein. The 

velocity potential is determined in the form of a double 
integr^, by application of the Fourier transform to the., 
boundary-value problem for this potential, and expressions 
for forces and moments are given in 'terms of this double . . . 
integral. The reduction of the int^ai expressions of 
reference 3 to forms desirable for flutter calculations — that . _ 
is. chordwise section forces and moments — is not given. 

SYMBOLS 

<t> disturbance-velocity potential 

x,y,z rectangular coordinates attached to wing 

moving in negative r-direction 
t ,57 rectangular coordinates used to represent 

space location of sources or doublets in 
ry-plane 

Z„ function defining mean ordinates of any 

chordwise section of wing such as y=yi as 
shown in figure 1 

w{x.yi,t) vertical velocity at surface of wing along 

chordwise section at y=yi 

Xo abscissa of axis of rotation of wing (elastic 

axis) as shown in figure 1 ^ 

t time 

h vertical displacement of axis of rotation 

Ao amplitude of vertical displacement of axis 

of rotation, positive downward 

armonlcally OsdDacInE Thin Rectangular Wings In Supersonic Potential Flow" by Charles 


E. Watkins, 19S0. 
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a 

angle of attack 


amplitude of angular displacement about 
axis of rotation, positive leading edge up 

h,a 

time derivative of A and a, respectively 

V 

velocity of main stream 

c 

velocity of sound 

M 

free-stream Mach number (F/c) 

is=Vm=*-i 



functions defined with equation (7) 


function used to represent space variation 
of source and doublet strengths 


function used to represent time variation of 
source and doublet strengths 

(a 

frequency of oscillation 

- APu 





k 

reduced frequency (coJ/F) 

i2=/3V(’? — 

i—v) 


represents functions of u, x, and M, defined 
in equations (15) 


represent functions of x, Zo, and 3, defined 
in equations (19) 

D. 

function used to denote doublet distribu- 
tions (see equation (22)) 

F. . . 

function defined in equation (28) 

G. 

function defined in equation (29) 

p 

density 

Ap 

local pressure difference measmed positive 
downward, defined in equation (31) 

b 

half-chord 

s 

half-span 

A 

aspect ratio (s/b) 

F 

total force acting on wing defined in equa- 
tion (32) 

Xi,is,XajXI 

components of total force coefficients, 
defined in equations (35) 

in 

total moment acting on \ving, defined in 
equation (36) 


components of total moment coefficients, 
defined m equations (38) 

P 

section force (total force at any spanvise 
station), defined in equation (39) 

Li,Lt,Li,Li 

components of section force coefficients, 
defined in equations (41) and (42) 


section moment (total moment at any span- 
wise station), defined in equation (40) 

iV/l J -214"4 

components of section moment coefficients, 
defined in equations. (43) and (44) 

Fn,6n 

functions related to and <?*, defined in 

appendix 


ANALYSIS 

BOUNDARY- VALUE PROBLEMS FOR VELOOTY POTENTIALS 

Consider a t.liin flat rectangular wing moving at a constant 
supersonic speed in a chordwise direction normal to its lead- 
ing e^e as shown in figure 1 . The boundary-value problems 
for the velocity potential for such a wing may bo conven- 
iently classified into two types associated with the nature 
of the flow over different portions of the wing. On the 
portion of the wing between the Mach cones emanating from 
the foremost point of each tip (region N in fig. 1 (a)) there 
is no interaction between the flow on the upper and lower 
surfaces of the wing. The type of boundary-value problem 

/ 


(a) 




FiovRE 1.— Sketch lUoatiatlng cboeen ooordlnale sTStem and the tvo dcfreofl of freedom 

a and h. 
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for this portion of the ^\Tng is referred to herein as “purely 
supersonic” and the velocity potential for region N is de- 
noted by <j>y. On portions of th.e wing within the Mach, 
cones emanating from the foremost point of each tip (r^ons 
Ti, Tj, and Ta in fig. 1 (a)), there is interaction between the 
flow on th.e upper and lower surfaces of th.e wing. The type 
of boundary-value problem for these portions of the wing is 
referred to as “mixed supersonic” and the velocity potentials 
for these regions are designated by and 0r,i respec- 

tively. The complete velocity potential at a point may then 
be expressed as 4>y, or 4>t^ according to th.e region 

that contains th.e point. 

As customary in linear theory, as applied to thin flat sur- 
faces, the boundary conditions are to be ultimately satisfied 
by the velocity potentials at the projection of the wing onto 
a plane (the ij^-plane) with respect to which all deflections 
are considered small and which lies parallel to the free- 
stream direction. Thickness effects are not taken into 
accoimt; hence, the velocity potentials are associated only 
with conditions that yield lift and are consequently anti- 
symmetrical with respect to the plane of the projected wing. 
It is therefore necessary to consider the potentials at only 
one surface, upper or lower, of the projected wing. The 
upper surface is chosen for this analysis. 

The differential equation for the propagation of small 
disturbances that must be satisfied by the velocity potentials 
is (when referred to a rectangular coordinate system x,y,z 
with the sry-plane coincident with the reference plane and 
moving uniformly in the n^ative r-direction, fig. 1) 


1^ 

c* 





( 1 ) 


The boundary conditions that must be satisfied by the veloc- 
ity potential are: (a) In regions Ti, Tz, Ta, and K the flow 
must be tangent to the surface of the wing or 



=w{x,yi,t)= 




( 2 ) 


For the particular case of a wing independently per- 
forming small sinusoidal torsional oscillations of amplitude 
lool and frequency os about some spanwise axis Xq and small 
sinusoidal vertical translations of amplitude \ho\ and fre- 
quency os, the equation of Za is 

Z«=e*“‘[ao(x— a:o)+Ao]=o<a:— aioD+A' (3) 

Substituting this e:q)ression for Zm into equation (2) gives 

w{x,yi,t)= Va-{-6ix—x^ +h (4) 

The velocity potential may thus be expressed as the sum 
of separate effects due to position and motion of the wing 
^sociated with the individual terms in equation (4) as 

( 6 ) 

DERIVATION OF 0.T 

The boundary-value problem in the purely supersonic 
region (fig. 2 (a)) is the same as that for the two-dimensional 
wing treated in reference 4. This problem is there shown to 
be satisfied by a distribution of soui'ces referred to, in this 
case, as moving sources because of the uniform motion; 
that is, 

(#).v(x,y,2,f )=— dv (6) 

In equation (6), TT'($,ij) represents the space variation of 
source strength and must be evaluated in accordance with 
the individual teuns of equation (4), and <pi is the potential 
of a moving source situated at the point (^,i?,0) that may be 
expressed as 

where w(t) is the time variation of source strength and the 
symbols with subscripts appearing in equation (7) are 
defined as 

jSc 


where Zm is the vertical displacement of the ordinates of the 
surface of any chordwise section of the wing (see fig. 1 (b)). 
(b) In regions Ti and Ti the pressure must fall to zero along 
the wing tips and remain zero in the portion of the Mach 
cones emanating from the foremost points of the wing tips 
not occupied by the wing. (Another condition, that the 
potential must be zero ahead of the wing and in the region 
off the wing adjacent to the Alach cones emanating from 
the foremost points of the tips, is automatically satisfied 
by the type of source and doublet synthesis employed in 
the solutions.) 


d/(r — g) , ■s/jy—vOirii—v) 

Cj3" j3c 

ni=y+^ V(z— S)*— 13*2* 

The time variation of source strength ia(t) for harmonic 
oscillations may be written as 


( 8 ) 
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(a) Purely sapersonlo region. 

(b) Mixed supersonic region. 

Fiovm 2.— Bketeh Illustrating areas ot Integration lor purely suixrsonlc and mixed auper- 
s<Hiloregl(Bi3olflon'. 

The numerator in equation (7) thus becomes 

wit— rO+mCt— 4-gf«(i-rs) 

—2e*'‘‘e ” * cos w— (9) 

Substituting equations (7) and (9) into equation (6) yields 


glut pX- 

4>w{x,y,z,t)= 

Jo 


\ ^ 

J.1 if 




(10) 

where, for briefness, 

_ uM ^Pu> 


and 




R= /(a: — —fiKy— i?)*— jS \(v—Vo(v 2 —v) 

The values of TT(fj’)) associated with the different terms of 
equation (4) are 


For A 

For Va 
For a(a:— Jio) 


Tr(i,»/)= 


W{^,r,)= Fao 


Ao 


a’o) 


( 11 ) 

( 12 ) 

(13) 


If any of the values of W{^, rj) given in equations (11), (12), 
and (13) is put into equation (10), the integration with 


respect to 17 can be readUy performed and the remaining 
integral evaluated as a series of Bessel functions. (Sec, 
for example, reference 4.) However, in order to be con- 
sistent with and to lead naturally to a succeeding part of 
the analysis the integrand is expanded into a Maclaurin’s 
series with respect to u. The expansion yickls 

giui rx-Pz (‘nx r/ 1 P 

<t>!f(x,y,z,t)=- — ir(i,n)[^(^aoi^ + aii^+ • • • 

+ • ■ -^ + (fflojf?+Ol2f^+ • . • 

+ a.2f»-ff+ ...)+•• .+(ao«f?*’"-'+a, 

-f. . . ...)+••■]«(’! «(^ (14) 

where tlie coefficients a„m are functions of w, *, and .1/; 
those coefficients involving w, up to and including the third 
power, are 

ffloi = 1 — ^ w a: - y rH J-* 
fw* , 

an=tw+« a; — ^3" 


W . tfW 

021=— y+-^ar 


081= ~T 


to>° 


iw’ 


Ou = 


2M* 


(16) 


Observe the following identity that is valid regardless of tlic 
highest power of w considered and that will be of subsequent 
use, namely 

Ooi+a;au-{- . . . -fa:’‘cf,iBl (16) 

It wiU be noted in equation (14) that the potential of a mov- 
ing source when expanded in terms of the frequency appears 
as a senes of terms similar to steady-state source potentials 
plus series of terms involving various powers of R. By 
grouping the terms in equation (14) ^vith respect to powei's 
of f, the. following form of the source potential convenient 
for later use is obtained : 

ptat px-Pt pV 2 r/ 1 

<l>N{x,y,z,t)= — ir(?,7?) |^(^aoijj-|-ao 2 ^+ • . . 

+ ao«-S*’“~®+ . ■ •^+J^aii^+ai2S+ • ■ • 


. -hr (a.ij+a. 
-h. . . .'jdnd^ 


.if 


(17) 


With the terms of the series grouped in tliis manner, in view 
of the fact that the differential equation (1) is independent 
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of I, it is apparent that the coefficient of each power of f in 
equation (17) is a solution to the differential equation. 

If the values of in equations (11), (12), and (13) 

are put into either equation (14) or equation (17), the inte- 
grations of each term can he easily carried out in closed 
form. Moreover it can readily be shown that, when all the 
terms involving S up to a given power are taken into account, 
the differential equation (1) is satisfied to the highest power 
of w considered. The boundary condition of tangential 
flow as expressed in equation (4) is satisfled exactly and 
does not depend on the order of « considered. 

Putting the values of "Pr( 5 ,i 7 ) in equations (11), (12), and 
(13) successively into either equation (14) or equation (17), 
carrying out the indicated integration, and setting 2=0 
yields for the velocity potential, to the third power of S at 
the upper surface of the wing, in the purely supersonic 
region: 

(18) 

where 


/,=:r-f.*-^,(2^*+3):r*+^(2^*+5)x^ 
/*=|(z-2x,)-^V-3a:o)-^*(2/S^+3)r»(a;-4xo)J 


Zca“ 


y ( 19 ) 


DERIVATION OP AND 


In order to satisfy the boimdary-value problem in regions 
of mixed supersonic flow it is convenient to start with the 
potential of a moving doublet. Then, for a given order of 
the frequency of oscillation, this potential, as will be shown 
in the following analysis, can he modified so that when 
integrated over the appropriate region the results will satisfy, 
as in the purely supersonic case, the differential equation to 
the given order of the frequency and wfil satisfy the condi-. 
tion of tangential fiow exactly. The potential of the type of 
doublet required may be obtained by partial differentiation 
of the potential of amoviogsource (see integrand of equation 
(17)) with respect to the direction normal to the plane of the 
wing, namely 

^ ^ + ^02^+ • • • + 

( 20 ) 

Examination of equation (20), like equation (17), shows 
that the coefficient of each power of f satisfies the differential 
equation and, since the differential equation is linear, it is 
permissible in synthesizing the solution to tiie boundary- 


fAlthoagb. the derlTatloa of these potentials In. NAOA TN 2064 led to correct results, tl» 
procedure ft^wed therein Is based on enoneous argommts. The present procedure Is 
correct and general. 


value problem to weight these coefficients separately. 
Furthermore the coefficient of each power of f consists of a 
term that has the form of a steady-state doublet potential, 
namely 

plus a series of other terms involving various powers of R. 
In the following analysis, attention is first directed to the 
treatment of the first term of the coefficient f", (expression 
(21)). The other terms wfll be treated subsequentiy. 

Expression (21) has the form that in steady flow is con- 
venient for treating the (antisymmetric) problem of satis- 
fying the condition of tangential flow for a distribution of 
normal velocity prescribed, at the wing surface, independent 
of 2 / but proportional to x“; that is, a weight or distribution 
fimction Z)„(f,ij) can be determined so that 

J £>«(^,n)^ n=(0,l,2, . . .) (22) 

where the region of integration r is the portion of the wing 
situated in the fore cone emanating from the field point 
(x,y,z) (shown in fig. 2(b) for the rectangular wing with 
2 = 0 ). 

The distribution function for rectangular wings may be 
easily determined when Da, the distribution function for this 
wing at constant angle of attack in steady supersonic flow, is 
known. The expression for Da is derived in reference 1 and 
found to be 

[V/3v (5-j8v)+ f sin-i-y/^] (23) 

From this expression and equation (22) it follows by direct 
substitution and reduction that 

Dn=nlf f-'-f Da{k,ii){dQ’'=n{ (f— X)’*"‘Do(X,ij)d\ 
Jo Jo Jo Jo 

(24) 

"With Dn known so that equation (22) is satisfied it may be 
demonstrated, with use of identity (16), that these “doublet- 
type” terms alone satisfy the condition of tangential flow. 
For example, let it be required to satisfy this condition for 
vertical translations; then, 

=h{aoi+xan+ . . . -)-x”a«0=^ 

The next step in the analysis is to consider terms of the type 

m>l (25) 

appearing in the coefficients of (equation (20)). It is to be 
noted that, when the power of w to which the potential is to he 
derived is chosen, the number of terms in equation (20) that 
are to be treated is determined by the expressions that 
contain u to the chosen order (see equations (15) for S to the 
third power) . 
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If the distribution function Z?„ (equation (24)) is intro- 
duced into equation (20) and the resulting equation is inte- 
grated over the appropriate region of the rectangle (fig. 2(b)), 
it is found that in the limit as z-^0 teims of the type 
given in expression (25) do not contribute to the resulting 
potential, but they do contribute to the resulting vertical 
velocity. Therefore, since it has been shown that the doublet- 
type terms, taken alone, satisfy the condition of tangential 
flow, the distiibution of vertical velocity would now contain 
extraneous terms that need to be canceled. This canceling 
may be achieved by what is essentially an iterative process: 
the functions! I?„ (equation (24)) are used with terms in 
equation (20) that have the form given in expression (26) to 
calculate the velocity that is to be canceled ; then, terms in- 
volving 1/R that may be individually weighted and which, 
of course, must satisfy the differential equation to the order 
of 3 to which the velocity potential is to be derived are added 
to equation (20); finally, by making use of equation (22), 
weight or distribution functions for these additional terms 
may be determined so as to cancel the extraneous velocity. 

This process is illustrated for <o to the third power as 
follows: 

where it k noted in the integrand that two extra terms have 
been added to equation (20) to accomplish the desired 
canceling. Substituting the values of TJT?i>7) defined m 
equations (11), (12), and (13) mto equation (26) gives the 
expression for the velocity potential at the upper surface of 
the wing 

r«[2F.-(2«+^..-^ *•)?,- 

j^24-2i3(a:— a:o)+ 2^8 (a:*— 2 a:zo)J F 3 + 

(x-Xo)j (2^=+ 1)F, j) (27) 

where the terms are grouped conveniently bj’^ the definition 
of Fn in the following integral: 

F„= sin~^ dx (n=l,2,3,4) (28) 


CThe functions F^, given in equation (28), and certain 
related functions are of particular importance in the remain- 
der of this development. Integrated values of this function 
for the first few values of n and expressions for related 
functions needed later in this anaij’sis are given in the 
appendix.) 

Examination of equation (27) shows that along the Mach 
line x=fiPf separating region Ti from region N, tho o.xprcssion 
4>t^ reduces to the expression for <f>if given in equation (18). 

The "corresponding potentials for regions Tj and T* can 
now be obtained. The potential is obtamed from equa- 
tion (27) by merely substituting 2a— y for y in equation (28) 
so that 

6n—j^ a:"~* sin~^ \/ ~" g ('i=l>2,3, 4) (29) 

The potential in region T* (that is for l^AfiK2) is a simple 
superposition of the potentials for regions N, Ti, and Tj, as 
in the steady case (see, for example, reference 5), and may be 
written as 

(30) 

FORCES AKD MO.MENTS 

Two t 3 TJes of force and moment coefficients are derived. 
First, in order to gain some insight into the over-all effect of 
aspect ratio on the forces and moments, expressions for total 
force and moment coefficients are derived. Tlien, in order 
to present expressions that are more suitable for use in 
flutter calculations, e.xpressions for section force and moment 
coefficients for any station along tho span are derived. 

Total forces and moments. — ^The local pressure difference 
between the upper and lower surfaces on the wing may be 
witten 

ij.= -2p(H+F||) (31) 

In order to derive expressions for total forces and total 
moments, it is only necessary to consider the velocity 
potential in two regions — either regions N and Ti or regions 
N and Tj. Therefore the expression for the total force, 
positive downward, on the wing may be writ ten as 

P=^2 j^jAptfdydx—2j^ jApr^dydx (32) 

where Apa is to be calculated from equation (18) and the 
integrations in the first term are to be extended over the 
shaded portion of region N shown in figure 2 (a), and wliere 
Apr, is to be calculated from equation (27) and the integra- 
tions in the second term are to be extended over region Tj. 
(The integrations in the first term are simple and may be 
performed by inspection. Those in the second term may be 
readily performed by making use of the relations given in 
the appendix.) 
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After the indicated integratLons have been performed and 
ail position coordinates involved have been referred to ^e 
chord 2b (but the original coordinate symbols maintained), 
the results can be written as 

{Li+iL^+a,{rz+iL^)'j (33) 

where the reduced frequency k is related to a and a by the 
relations 


j b(j> 6j3^ _ 


and where 


^-W'-w {si*- 

^ [(3/3*+4)+4i3^2a(2+i3*)l j 


2:^=^03*-l-2^^Xo)+^(5+r+12^^ro)+ 

(2 + 3M-^ (8+4^*+20/3*a:o+5^^£o)] 

(35d) 

The quantities Li (i=l, 2, 3, 4) are the in-phase and out-of- 
phase components of the total force coefficients, Li and £? 
being the in-phase and ij and i* being the out-of-phase 
components. It will be noted that Li and Li are associated 
only with vertical translations of the wing and are inde- 
pendent of axis-of-rotation location Xo- The components 
Zj and Li are associated with angular position and rotation 
of the wing about any axis x—xo and depend partly on the 
location of Xo- 

The total moment, positive clockwise, on the wing about 
the arbitrary axis of rotation x=Xo is 

ai= ~2 J'^J' (x—Xo)A.pjidydx—2j'^ J (x—Xo)Apr^dp dx 


If steps similar to those required to obtain equation (33) 
are performed, there is obtained 

]S:=--8p 

(37) 

where 

3Ii=^(2-3xo)-^J3|3*-l- 1 -4xo(2-b)3^] (38a) 


“ 2 ro) - ^ [3 + iS"+ 4xo(j3*- 1) - 8/3*Xo"l - 
^ { (2 - 3 Xo) - ^J4 (4 -b 3 /3^ + 5 Xa(3 |S*+ 3 18»- 4) - 
20|8V(|S*+2)]| (38c) 

321=^ [203*- 1) -3xo(2/3^- 1) -h 6j8’-Xo=] -h 

(20 + 4j8*-2 5xo-l-40^*Xo- 60/3*Xo*)-F 



[2 0(2 + j3=) -24xo(2 -3^*-,3^ -30|S*Xo*(4-[-i3*)l I 

(38d) 

The quantities Mi and 3^ are, respectively, the in-phase and 
out-of-phase components of total moment coefficients about 
the axis x=Xo associated with vertical translations of the 
wing; Ms and Af* are the corresponding components due to 
angular position and rotation of the wing about x=Xo. 

It is of interest to note in equations (35) and (38) that the 
components Li and Mi do not involve the reduced frequency 
k. The effect of frequency on these two components comes 
from terms involving the frequency to the fourth and higher 
powers; but for values of k thought likely to be encountered 
in supersonic flutter (i<(0.1), the contribution of these higher- 
power terms to any of the components in equations (35) and 
(38) is, for the most part, negligible. 

Section forces and moments. — ^The section forces and 
moments at any spanwise station are derived by integrating 
the pressure difference along the chord for the forces and the 
pressure difference multiplied by a moment arm for the 
moments. Since the distribution over the entire wing is 
symmetrical with regard to the midspan section, it is only 
necessary to derive expressions for the forces and moments 
at any station of the half-span adjacent to the origin. (See 
figs. 1, 2, and 3.) 

Under the restrictions previously stipulated, two cases 
that can arise are considered (see fig. 3) : (1) the Mach lines 
from the tips do not intersect on the wing (or M^>2), and 
(2) the Mach lines intersect on the wing but the Mach line 
from one tip does not intersect the opposite tip ahead of the 
trailing edge (or lgAj8^2). Only the final forms of the 
section force and moment equations are given. These 
forms are easily calculated by deriving the pressure difference 
for the different regions from the appropriate velocity 
potential, making use of figmre 3 to determine the limits of 
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integration for the regions involved, and using the relations 
given in the appendix to carry out the more troublesome 
integrations. The integrated expression for any region can 
then be reduced to the forms 


P=-4p6FW-' [^(ii+iXs)+oo(i3+ii4)] (39) 

and 

M„ = -4p6*FW“' (40) 


where the position coordinates are referred to the chord 
length 26. 

The components of force and moment coefficients for the 
half-span adjacent to the origin are as follows: 

Case 1 (see fig. 3(a)): For any section between the tip 
and the point where the Mach line intersects the trailing 


edge, or where the components of section force 

coefficients are 




=-A( 




Pi 


) 


r 4 r 1 = 5 6^^+3|S’-l = 

^ P,-(1.^2ro)Pi+ P»- 

. /3* 


(41a) 

(41b) 

(41c) 


2^%-4|3^a:o)F,-6^*(4-3ro)K-(l-7^'-20/3^)F4] | (41d) 


where P» (n=l, 2, 3, 4), given in the appendix, is obtamed 
from Fn (equation (28)) when x=2b. For any section 
between the point where the Mach line intersects the traU- 


1 A 

ing edge and the midspan, or where —^ 3 / the compo- 

fS £i 


nents of section forces are: 



^ 1 JIP* 

W [(3+/3*) + 6/3*Xo] 


r (42) 


1 

^4=;^ [(|9*- 1) - 2/3*Xol + '^(5 -t- 2/3%) J 




(a) Mach lines from tips do not Intersect on whiB. 

(b) Mach lines from tips Intersect on wine but Mach line from one tip does not Intersect 

opposite tip. 

(c) Mach lines bom tips Intersect on wing and Mach lino from one tip Intorsi’cls opposite 

tip at traQing edge. 

Fiocbe 3.— Sketch Illustrating different Mach line locations accounted for In onali'Sls. 

The components of section moment coefficients for case 1 
are as follows: 


For 0 < 3 /< 


(1-2xo)F+2xo 


ii!±lp 3<3H2 

iS* i3* 



(43a) 


As a check on the results in equations (41) and (42) the 
expressions in equations (41) reduce to those in equations 

(42) when y=i- 




4j3»+l 



(43b) 
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M,=|;;[p 3^o+3xo^ 

(6/9^+3/3*-l)(l-2io)+4/3W(l + 2^‘Oi , 
p ^s+ 

6(l+|3^_s 20i9^+21^*+3 5l 

p ajoi'a 3^4 ■t'^J 


(43c) 


+ 


^^{^[(4/3^-2i3Hl)(l-3xo+3cro*)+l-3jro^F2- 
(2-4j:o+3jo’)^,-^ 4 (8i8^+4i3^-l)f4+ 


j^(S|3’+3)#.}) 


(43d) 


aadfor^gi/gy. 


^(3-4aro) 


— 3^3 (2 3io) 

(1 — 2xt^~ 

M,=^, (1-2 jo)-^3 (3+(3*-4ro+4/3^a:o-8|3^Xa^ 
[2(^^-1)-3 xo(2^*-1)+6M + 


APk 

150 


^ [4(5+^^+5ao(8^^-5)-60/3V] 


y (44) 


The expressions in equations (43) reduce to those in equa- 
tions (44) when «=i. The expressions in equations (42) and 
P 

(44) correspond to the more exact two-dimensional compo- 
nents of force and moment coefficients derived in reference 4. 
For values of i*<C0.1 these expressions yield, for the most 
part, values that are in good agreement with those that 
may be calculated from the tables in reference 4. 

Case 2 (see fig. 3 (b)): For any section between the tip 

at y=0 and the point where the Mach line from the tip at 

y=2s intersects the trailing edge ^or where 0 ^^7 

the components of section force coefficients are given bj' 
equations (41) and the components of section moment co- 
efficients, by equation (43). For any section between the 
point where the Mach line from the tip at y=2s intersects 

the trailing edge and the midspan ^ or where A—-<Cy^-^ji 

the components of section force coefficients are 

ih+SS\+^ (4oa) 

5:)+^[(2/3*-1)(#.+ (P3)-3|S* (#3-1- #a)]}- 




Bk 

21363T — 08 


^Pk\ 
/3® / 


(45b) 


(^i+54)-2(1-23:„)(F4+(?0-1- 

M!±|^ 

[(3-l-i3»)-h6«| (45c) 

i4=^{^[(2-xo)(#i+(50-^^ (#*+<^2)]-b 


APk 
3 jS® 


[3(1 -)S^-i-2i3^-f2^^aro-4^^j:a) (F3-f-<93)- 


6)SX4-3xo)(^’3-f-G'3)-(l-7i8*-20|8^)(#4-fe4)]|- 

{^[03^-l)-2|3^XaI+^(5-h)3*-M2i3=3;o)| (45d) 

where On (n=l, 2, 3, 4), given in the appendix, is obtained 
from On (equation (29)) when x=2b. The corresponding 
components of section moment coefficients are 

My j^[(l -2x,){fy+dy)A-'^^^^P^ (^2+^3)]- 


2jCo(2^^-H) 

(#3+5,)+^^ (2-3ro)} 




(46a) 




(l_ 2 r„)_:^ (3-4xo)] 

(46b) 

(6j3^-h3i3^-l)(l -2xo)-h4|3^Xo^l -f 2j3^ [F,+U,) + 

{^(l-2xo)-^ K3-t-i3^-4xo(l-i8^-8^*Xo^| (46c) 

M,=f^(j[{l-Xof{Fy+^y) Xoif,+^,] - 

^^(f3-h^j}f^*{^4[(4^^-2^=+l)(l-3a:o+3xo^ 
1 - 3xo^ -(2 -4 xo+3xo=0(^3+^3) - 

(8|3*-f 4^=>-l)(f4-t-54) (5|S*-h3)(f3-h^J I)- 

{ 3 ^ [2(/3^-l)-3xo(2i3*-l)-l-6j3^Xo^-h 

^ [4(5+^^-f-5xo(8^*-5)-60|3^Xo^| (46d) 


-48 
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(a) Li. 

Figure 4j-Components of total force ooefflclents aa functions of M>foru-0.4, t-0.02, and 
TRiious vEdaes of A. 



Figure 4.— Continued. 


For the limiting condition of case 2 — that is, when the 
Mach line from one tip intersects the opposite tip at the 

trailing edge, or -4=^ (sec fig. 3(c)) — the> components of 

section force coefficients are given by equations (45) and 
the corresponding components of moment coefficients, by 
equations (46) . 

SOME PARTICULAR CALCULATIONS AND DISCUSSIONS 


From the expressions for total force and moment . co- 
efficients (equations (36) and (38), respectively) the over-all 
effect of aspect ratio on the magnitude of the forces and 
moments can be calculated for particular values of the 
parameters M, k, x^, and A. Examination of these equa- 
tions shows that varying some of the parameters might 
cause some terms in the equations to vanish and to change 
sign. For example, if *o is continuously increased from some 
value less than 1/2 to some value gi-eater than 1/2, thefii-st 


terms in the expressions for and vanish at a-o=^ and 


change sign when aj© becomes greater than 1/2. In particular, 
decreasing the aspect ratio decreases the components of force 
and moment coefficients Lj, ij, ilfi, Mt, and Ml but 
increases the two important components Li and Mi. 

Although the effect of aspect ratio may change considerably 
with only a small change in any one (or more) of the param- 
eters M, k, and x^, some insight into what the over-dl 
effect might be can be gained from calculations of all the 
components of total force and moment coefficients for various 
values of M and A and fixed values of the parametei-s t and 
aio- Results of such a set of calculations are presented in 
figures 4 to 7. 

In fiiguros 4 and 5 the components of total force and 
moment coefficients for various values of A and for ro=0.4 
and A:=0.02 are plotted ^ functions of M*. The curves in 
these figures calculated for infinite aspect ratio correspond 
to the two-dimensional results of reference 4. Tho dashed 
curves represent calculations for aspect ratio and Slach 
number, combinations that cause the Mach lines from one 
tip to intersect the opposite tip at the trailing edge so that 
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4 


Cc) U. 

FicuaE 4.— Contlnned. 

along the dashed curves the aspect ratio is not constant but 
varies with according to the previously given expression 

■ 4 _ 1 _ 1 

/3 VilP-1 

The difference, at any value of iV/*, between the dashed 
curves and the curves corresponding to infinite aspect ratio 
in figures 4 and 5 is, therefore, for the chosen values of 
k and Xo, the maxiimmi effect of aspect ratio on the com- 
ponents of total force and moment coefficients for a rec- 
tangular wing under the restrictions of the foregoing analysis. 
It win be noted in figures 4 and 5 that, when the aspect 
ratio is small, the deviation of the three-dimensional results 
from two-dimensional results may be quite large. 

In figures 6 and 7 the components of the total force and 
moment coefficients are plotted as functions of aspect ratio 
for ro=0.4, k=0.02, and some particular values of M. It 
will be noted in these figures that all the components of 
force and moment coefficients undergo rapid changes with 



Figuke 4.— Concloded. 


respect to varying aspect ratio when A becomes less than 
4 or 5. It may be remarked that the directions of the changes 
with respect to aspect ratio appear to be such that they would 
have favorable effects on the flutter characteristics of a wing. 

The spanwise distribution of the components of section 
force and moment coefficients computed from equations (41) 
to (44) for A=4, Xq=0A, I’=0.02, and M=2 are plotted 
in figures 8 and 9. The portions of the curves in these 


figures corresponding to values of in therange — ~ 


are the two-dimensional values, and the effect of aspect ratio 


may be noted in the tip regions, Ogy^^and 

as deviations from these two-dimensional values. 

In conclusion, it may be stated that, in regard to the 
effect of aspect ratio on supersonic flutter, an important item 
that has not been discussed herein but can be studied for any 
particular case with the aid of equations (33) and (36) 
is the change in center of pressure, associated with prescribed 
motions of the wing, with change in aspect ratio. An 
investigation to find the effect that thickness might have 
on the center-of-pressure location is also needed. An 
extension of the foregoing analysis to include the effect of an 
aileron as an additional degree of freedom would follow 
in a straightforward manner. 


Langlet Aeronautical Laboeatort, 

National Advisory CoMinTTEE fob Aeronautics, 
Langley Field, Ya., January 6, 1960. 


























APPENDED 


SOME INTEGRATED VALUES OF AND OTHER RELATED FUNCTIONS 

Values of F„ and Q „, — The values of the functions F„ (equation (^)) and (?„ (equation (29)) for the first few values 
of n are as follows: 


a:" ‘sin“‘ -y/^ylxdx (71=1,2,5,...) 
Fi=^fiy(z—^y)+xsm-^ 

Fi= ^'^^'^ ^/Py(x — <3y)+y sin“‘ -y/py/x 

c, 3a:’*+4j3j/a;+8/3V nr-, . a:’ • ‘-tt-t 

^3= Jg ^ V^ 2 /(a;— /3i/)+-g sm 'vi3y/a: 

Ft= ^ ^ V/3y(a:— (32/)+-^9m ^^J^yfx 

sin- 

(n=l,2, 3, . . . ) 




'K^s-y) 


(?) = Vl3(2s— y) [a; - i9(2s — y)] + a: sin-‘ ^ 

g_ a:+2ff(2ar-y) ^^^(^^-y)[^_^(2^_^)]^|sm-^ ^^^ — - 



+ ■,l ^i2s-y)[x-K2s-y)] sin- - /IB 

1 57 6 5 V ^ 

Values of F^ and — The following expressions define the functions and §n appearing in equations (41), (43), 
(46), and (46) in the body of the report. In these expressions the variable y has been referred to the chord 2b; that is 

y/2i has been replaced by y and in the expressions for the ratio s/6 has been replaced by A: 

^i = Vi8j/(l— j8j/)+sin“S/i^ 

V<3y(r— ffi/)+|sin-^ 

^ yPy(l ~^y)+| sin-^ 

+40^,+48yv-y 4^ Y+ 1 28^^. ^ 
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(?i=Vl3(A— y)[l— j9(^— y)]+sin.-Vi3(^— 3/) 
§_1±2P^:^ -fi(A-y)] +1 sin-^ M^) 




45 

g__ 35+40«^-,)+48^-(^-V)’+64m-tf+12Sm-rt‘ V ^(A-y)[l-W^-i,)l +l4m- 

15/5 5 

Some integral relations for F„ and (?,. — Some pertinent integral relations for are as follows: 


F^dx=xF^—F„ 


xF„dx=-^(x^Fn—F^ + 2 ) 


X 
X 

j^^x^F,dx=^(x^F,-F„l,) 

{x—x,^F„ 

£ (x - XoYF, ^ ^ i^»-a:o*^;+i+a;o>,+2-|- 


dx=^^-^ F„+XoF„+,-^ F„+, 


Corresponding integral relations for (?a may be obtained from these relations by simply replacing J* by <?. 

Integral relations for F „, — ^Integral relations for F„ that may be used in calculating total forces and moments from 
sectional forces and moments are as follows: 
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